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CHARACTERISTIC  EMISSIONS  FROM  AN  UNDERGROUND  EXPLOSION 
Second  Seml-Anniial  Report 


I.  INTRODUCTION 

The  investigation  of  characteristic  emissions  from  an  underground  nuclear 
explosion  described  in  this  report  represents  a  portion  of  the  total  effort 
initiedly  undertaken  by  the  staff  of  the  Engineering-Physics  Company  in  June  1961. 

A  technical  summary  of  the  work  carried  out  prior  to  this  reporting  period  is 
contained  in  the  first  semi-annual  report  issued  January  29,  1962. 

The  VELA-UNIFORM  program  was  initiated  for  the  purpose  of  developing 
techniques  and  equipment  for  the  detection,  at  a  distance,  of  undergrovmd 
(  nuclear  explosions .  One  of  the  more  difficult  problems  of  this  task  is  to 

\  effectively  distinguish  between  the  seismic  disturbances  caused  by  an  underground 

nuclear  explosion  and  those  caused  by  earthquakes.  For  the  purpose  of  developing 
'  .  both  equipment  and  techniques  for  detection  it  is  important  that  the  characteristic 

wave  shape  of  the  mechanical  disturbance  be  known. 

;  Under  sponsorship  of  the  Defense  Atomic  Support  Agency,  EPCO  has  been 

given  a  responsibility,  in  connection  with  the  VELA-UNIFORM  program,  to  determine 
the  shapes  of  seismic  waves  emitted  from  the  region  close  to  an  underground  nuclear 
explosion.  The  wave  shapes,  given  in  such  experimentally  measured  quantities  as 
particle  acceleration,  velocity,  and  displacement  and  expressed  in  terms  of  the 
properties  of  the  medium  surrounding  the  explosive  source  are  of  particular 
interest  to  those  having  the  responsibility  for  determining  whether  the  seismic 
signals  received  are  generated  by  a  nuclear  explosion  or  by  an  earthquake. 
Information  on  wave  shape  will  also  be  of  interest  to  those  designing  experiments 
involving  underground  explosions. 
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Tht  nature  of  tha  aaohanloal  waves  emitted  from  the  so  called 
non-linear  region  close  to  the  explosive  source  has  heen  viewed  In  this  report 
as  the  following:  the  bomb  Is  detonated  and  the  large  thermal  energy  release 
vaporises  the  bomb  parts  as  well  as  some  of  the  soil  In  the  iomedlate  neighbor¬ 
hood  of  the  source.  Because  of  the  disparity  between  the  properties  of  the 
gas  within  the  cavity  and  the  soil,  disturbances  caused  by  the  explosive 
source  reverberate  within  the  cavity  and  the  gases  soon  become  uniformly  dis¬ 
tributed  through-but.  The  exceedingly  large  pressures  within  the  cavity  cause 
the  cavity  walls  to  be  pushed  outward  rapidly  tending  to  form  a  spherical  room. 

As  the  walls  are  forced  outward  a  spherical  shock  wave  advances  ahead  crushing  the 
soil  as  It  propagates  radially  away  from  the  explosive  source.  The  crushed  soil 
can  be  thought  of  as  a  "spherical  piston"  which  Is  continually  gaining  mass  and 
absorbing  mechanical  energy  from  the  gas  In  the  cavity,  all  the  while  slowing  down. 

In  the  ease  where  the  shock  front  velocity  decreases  below  the  small  distur¬ 
bance  propagation  velocity  of  the  medium  a  new  wave  will  be  generated  and  this  wave 
will  propagate  radially  away  from  the  front  of  the  spherlcsd  shock  wave.  This 
disturbance  Is  sometimes  referred  to  as  a  sonic  precursor.  As  the  precursor  and 
the  main  disturbance  propagate  to  greater  distances  they  finally  reach  a  free  sur¬ 
face,  namely  the  surface  of  the  earth,  and  at  this  Instant  a  rarefaction  wave  begins 
to  propagate  back  Into  the  soil  tending  to  relieve  the  conpresslve  stress  of  the 
direct  wave  while  conserving  energy  by  particle  acceleration.  If  the  resultant 
stress  produced  by  the  combination  of  the  rarefaction  wave,  conpresslon  wave,  and 
the  llthostatlc  pressure  is  a  net  tensile  stress  In  excess  of  that  required  for 
rupture  of  the  soil,  the  earth  will  spall  and  one  or  more  slabs  of  soil  will  be 
hurled  Inward.  As  these  slabs  fall  back  due  to  gravity  they  Impact  with  the 
underlying  earth  generating  a  secondary  mechanical  disturbance  having  a  wave  shape 


INOINIIRINO-PHYSICS  COMPANY 


-  3  - 

very  different  from  the  direct  disturbance  Initieilly  generated  by  the  explosive  source. 
As  the  pressure  In  the  cavity  diminishes  the  cavity  wall  may  no  longer  be  able  to 
support  the  crushed  soil  at  its  exterior  and  the  cavity  will  collapse  setting  up 
an  additional  mechanical  disturbance  which  can  also  be  of  significant  magnitude. 

The  total  problem  with  which  we  are  concerned  can  be  divided  into  two 
not  altogether  separable  smaller  problems. 

1.  Generation  of  mechanical  disturbances; 

2.  Propagation  of  the  resultant  disturbances. 

An  investigation  of  the  generation  of  mechanical  disturbances  must  take  into 
account  not  only  the  form  of  the  initial  energy  of  the  source  and  the  way  in  which 
it  is  first  converted  into  mechanical  energy  but  also  conversions  from  one  form 
of  mechanical  energy  to  another.  For  example,  we  have  the  case  where  direct 
mechanical  waves  cause  a  moving  spall  which  in  turrv  upon  impact  gives  up  its 
energy  in  the  form  of  a  secondary  wave.  The  problem  of  propagation  of  mechanical 
disturbances  is  complex  when  considered  in  its  entirety  because  of  geometric 
considerations  such  as  lack  of  synsnetry,  irregular  surfaces,  layering,  and  so 
forth.  Additionally  the  various  types  of  geologic  material  and  their  non-linear 
mechanical  behavior  bring  in  the  requirements  for  rather  complex  constitutive 
equations.  The  resultant  eqiuations  obtained  when  analyzing  close-in  wave 
propagation  are  non-linear  and  of  a  form  which  will  admit  of  solutions  only  by 
numerical  methods.  A  complete  detailed  analysis  of  the  entire  wave  generation  and 
propagation  problem  would  involve  a  momumentea  computational  task  which  ^rauld 
more  than  exceed  the  practical  capacity  of  our  largest  present  day  electronic 
computers.  Hence  we  r-ust  be  satisfied  with  a  treatment  of  the  total  problem 


( 
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a  piecewise  fashion  making  realistic  simplifying  assumptions  when  necessary. 
Tliij  rnntiiiiits  nf  this  rsyn^  mi  in  f'll*  "  ~ 


y' 

tte  propagation  of  waves  very  close-in  to  the  source  of  the  erolosion  where  the 
soli  is  subject  to  very  high  compressive  force^  Results  of  the  analysis  are 
given  in  a  form  similar  to  typical  experimental  data — particle  velocity  and 
displacement  each  plotted  against  tlme^  shock  front  arrival  time,  and  peak 
particle  velocity  versus  radial  distance  from  the  explosive  source.  \  ts  ta¬ 
il  I  -III—  II  I  I  7  m  take 

into  aseount  uioi'B  deULTlHa  etj.iiaOieas  aff  stale  gild  Bll8fF^curate'''Sqttgtiens  of 

ghs  dlsetnislea  ef  QuLlluH'tll  lisale'^flth  the  transition  of  the  disturbance 


to  seismic  waves^  Results  are  expressed  in  terms  of  particle  acceleration, 
velocity  and  displacement  versus  time.  iiM0a«M«i»rBHHa-«amadidM»  ^he  elastic 
precursor  which  propagates  ahead  of  the  main  disturbance  .in  the  event  that 
this  main  disturbance  at  some  time  propagates  subsonically^yy  Some  qualitative 
results  of  the  analysis  are  presented  showing  the  wave  forms  expected  at  a 
hypothetical  gage  station.  The  manner  in  which  information  obtained  from  the 
precursor  can  aid  in  the  determination  of  equation  of  state  and  other  details 
of  the  wave  propagation  phenomena  is  also  discussed.  secondary 

wave  generation  as  the  result  of  spall  impact^  inwee**piaiiiBj^'^Mdefwater 
■pteeione  liaoe  long  been  aware  ef-the  'sgeufraaiy  wavtu  jiew^Satad  in  ■ooBB«otian.4*ith 
thr  nn  aallni  bulk  rnirt'ttiHisn  phenimijuj'iai  iirff*  iiiii  Iiit*  *^mtTr  fit 

B^ei'lments  ehow 

-I... 

®  f\ 

of  a  type  of  mechanical  disturbance  quite  unlike  the  usual  body  wave  disturbances 
issued  directly  from  the  explosive  source. 


asjik^hTinnin~n~n  1iT-TT~nnnlT  <r  •in'^Trcr — — 
a  rudimentary  analysis  which  indicates  the  generation 
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II.  WAVKB  IN  CRUBHABLE  MEDIA 

In  studying  vave  propagation  In  a  crushable  media  we  must  first  decide 
upon  an  equation  of  state.  Fig.  1  Illustrates  the  equation  of  state  to  be  used. 

We  assume  as  a  simple  first  model  that  the  medium  under  consideration  when  subject¬ 
ed  to  an  external  mechanical  force  crushes,  resulting  in  an  instantaneous  change 
from  the  Initial  density  to  some  final  density.  We  further  assume  that  this  final 
density  Is  Independent  of  the  applied  mechanical  force.  Dry  tuff  and  dry  alluvium 
behave  somewhat  In  this  fashion  and  while  the  assumption  that  crushing  occurs  In¬ 
stantaneously  Is  unrealistic  on  a  physical  basis,  we  will  find  that  this  assumption 
will  mainly  affect  the  shape  of  the  wave' s  leading  edge.  (This  assumption  Implies 
that  after  the  material  has  been  crushed  the  density  gradienx;  in  all  directions 
from  a  given  point  will  be  zero. )  We  begin  by  neglecting  the  energy  equation 
and  consider  only  conservation  of  mass  and  momentum,  thereby  neglecting  energy 
conversion  to  other  forms.  Going  back  to  the  equation  of  state  we  note  that  in 
Fig.  1  the  sloping  broken  line  represents  the  Jumn  in  pressure  and  density  across 
the  shock  front;  the  pressure  Jumps  from  to  P,  euid  hence  the  density  Jumps 
from  to  p.  Upon  relaxation  of  the  mechanical  force  from  P  back  down  to 
P^  the  density  remains  constant. 

We  consider  the  following  problem;  a  spherical  cavity  in  a  crushable  medi¬ 
um  of  infinite  extent.  We  neglect  the  wave  propagation  phenomena  in  the  expanding 
gases  and  other  materials  in  the  cavity.  We  assume  that  the  pressure  in  this  cav¬ 
ity  quickly  attains  a  state  where  the  pressure  is  uniform  throughout  the  cavity. 

As  a  force  exerted  by  this  pressure  on  the  cavity  wall  forces  the  wall  away  from 
the  center  of  the  cavity  we  assxime  the  following  relationship  between  pressure, 

P,  in  the  interior  gas  And  the  location  of  the  cavity  wall 
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dlaeaalonlesB  density^/p^ 

Fig.  1  --Solid  curve  Indicates  state  points  of  shocks  in  the  P-V  plane. 
Sloped  broken  line  indicates  shock  rise  fron  initial  pressure, 

P  to  peak  pressure^P  and  dloensionless  density,  p/p  .  Vertical 
broken  line  Indicates  fall  to  final  state  at  pressure, P^  and 
dimensionless  density,  p/p^. 
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(1) 


where  Is  the  Initial  pressure; 

Is  the  Initial  position  of  the  cavity  wall; 
a^  Is  the  position  of  the  cavity  wall  at  time,  t; 
y  Is  a  constant. 

Note  that  we  also  assume  this  wall  to  be  Impermeable  In  that  no  leakage 
of  the  gas  occurs  through  the  cavity  wall.  As  the  cavity  wall  is  forced  outward 
a  shock  wave  will  propagate  ahead  of  It.  This  shock  wave  will  be  in  the  form  of 
a  spherical  disturbance  such  that  on  one  side  the  medium  has  experienced  no  mech¬ 
anical  motion  and  that  at  the  front  the  material  receives  its  first  motion.  A- 
cross  this  front  wo  must  have  conservation  of  mass.  Taking  into  account  the  equa¬ 
tion  of  state  previously  described  we  have  for  the  conservation  of  mass  the  equa¬ 
tion 

(R^  -  r^)p^  =  (R^  -  ?^)p  (2) 

whore  R  is  the  position  of  the  shook  front  measured  from  the  cavity  center; 

r  is  the  undisturbed  position  of  a  particle  in  the  medium  behind  the  shock 
front; 

r  is  the  position  of  the  particle  which  was  initially  at  r; 
p^  is  the  initial  density; 
p  Is  the  final  density. 

Solving  Eq.  (2)  for  r  we  have 

?  =  [R^d  -  p,/p)  +r^(pyp)]^/^  (3) 


i 
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The  displacement,  u^,  Is  then  given  by 

=  r  -  r  (4) 

r 

e 

The  equation  for  the  particle  velocity,  r,  can  be  obtained  directly  from 
Eq«  (3)  by  a  differentiation  with  respect  to  t,  noting  that  the  only  function 
which  Is  dependent  on  t  In  Eq.  (3)  Is  R.  Differentiating  we  obtain 


?  =  ^  R(1  -  p>)  (5) 

In  Eq.  (1)  we  had  written  the  pressure  interior  to  the  cavity  in  terms  of  the  in¬ 
itial  pressure,  P^,  the  initial  position  of  the  cavity  wall,  a^,  and  the  position 
of  the  cavity  wall  at  some  time,  t,  a^. 

Using  Eq,  (1),  we  have  the  force  exerted  on  a  small  area  of  the  wall  of 
the  cavity  at  sane  time,  t,  given  by 


<P 


(6) 


where  cp  is  a  small  solid  angle  generated  from  the  center  of  the  cavity  and  hence 
the  impulse,  at  time,  t,  is  given  by 


I  =/%pP  aJ'S 
9  wf  ~  O  0  o 


3y-  (2  -  3y) 


dt 


(7) 


The  momentm,  transferred  to  the  material  which  is  moving  as  a  resvilt  of  the 
impulse  given  in  Eq,  (7)  is 


%-/ 


R  . 

rdm 


(8) 


where  dm  is  the  mass  of  a  segment  of  an  elemental  spherical  shell. 


INOINilRINO-PHYSICS  COMPANY 


-  9  - 


We  nay  Eq.  (8)  In  the  form 

^  (prpr^dr  (9) 

“o 

We  may  write  the  equation  for  momentum  in  a  more  useful  form  by  using  Eq.  (3)  end 
changing  the  variable  of  integration  to  obtain  the  expression 

H,  =  -  pyp,  (XO, 

o 

which  Is  valid  for  any  fixed  time,  t.  Looking  closely  at  Eq.  (10),  we  note  that 
this  Integral  can  be  simply  evaluated  and  we  obtain  the  final  form  for  the  momen¬ 
tum  given  by 

=  <i.py(x  -  pypjfp/p^  -  [(p/p„  - 1)  --(y  (XX) 

From  classical  mechanics  the  Impulse  must  be  equal  to  the  momentum,  hence 
we  may  equate  Eqs.  (7)  and  (11),  obtaining 

f»^(l  -  py»)  +  »o  (ep/pl) 

=  P>^(X  -  p^p)  ^p/b„  -  [^(p/p„  -  X)  +(^  J  ^^^)(X2) 

This  equation  is  a  non-linear  Integro-dlfferential  equation  but  it  can  be  solved 
quite  simply  by  nimerical  methods  upon  application  of  a  finite  difference  analogue. 
By  iterative  techniques  the  value  of  t  can  be  found  while  incrementing  R.  Using 
this  procedure,  values  of  R,  the  position  of  the  shock  front,  as  a  function  of 
time  can  be  determined.  Fran  this  result,  using  Eqs.  (3)  and  (4),  the  displace¬ 
ment  as  a  function  of  time,  and  from  Eq.  (3)  the  particle  velocity  as  a  function 
of  time  can  be  found.  A  subsequent  differentiation  of  Eq.  (3)  will  yield  an  ex- 
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preaaion  for  the  particle  acceleration  for  a  particle  originally  at  a  radial  dis¬ 
tance,  r. 

Ccaputer  programs  were  written  for  an  IBM  1620  cosymter  to  solve  Bq.  (12). 
The  Inputs  to  these  computations  were 

1.  Initial  density  of  uncrushed  soil,  p^. 

2.  Ratio  of  density  of  crushed  soil  to  density  of  uncruahed  soil,  p/p^. 

3.  Initial  cavity  radius,  a^. 

4.  The  gas  constant,  y. 

5.  The  pressure  In  the  cavity  at  time  zero. 

6.  Distance  from  explosive  source  to  five  hypothetical  gages  distributed 
along  a  radial  line  from  center. 

In  Fig.  2  we  have  Illustrated  the  cavity  wedl  and  the  five  hypothetical  velocity 
and/or  displacement  gages  at  three  different  times.  The  position  of  the  shock 
front  Is  also  indicated  In  each  case.  Fig.  3  shows  the  displacement  versus  time 
at  each  of  the  five  hypothetical  displacement  gages  for  the  data  described  in  the 
caption.  Fig.  4  shows  particle  velocity  versus  time  at  each  of  the  five  hypotheti¬ 
cal  velocity  gages  for  the  same  case.  Fig.  5  shows  the  shock  front  arrival  time 
for  a  crushable  soil  under  the  same  conditions.  Fig.  6  Is  a  plot  of  particle  ve¬ 
locity  versus  distance.  This  plot  Is  again  made  on  log-log  paper  so  as  to  clearly 
Illustrate  the  character  of  the  dependency  of  velocity  on  the  exponent  of  the  dis¬ 
tance,  r.  Eq.  (5)  shows  us  the  relationship  between  Figs.  5  and  6,  Indicating 
that  the  particle  velocity  Is  directly  dependent  upon  the  shock  velocity  and  posi¬ 
tion  of  the  shock  front.  From  further  calculations  It  appears  that  the  shock  front 
velocity  decrease  with  radial  distance  Is  affected  largely  by  the  assumed  value 
of  the  constant  exponent  of  the  gas  In  the  cavity  (see  Eq.  (1)).  Our  very  simple 
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Curves  A,  B,  C,  B,  and  E  represent  respectively  the  velocity  as  a  function 
of  time  for  each  of  the  five  hypothetical  gages. 
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model  here  aaeumee  that  the  density  of  the  crushed  material  Is  Independent  of  the 
pressure  level  at  the  shock  front;  while  for  crushable  materials  this  appears  to 
be  very  nearly  correct,  we  know  from  experimental  evidence  that  In  most  cases  the 
crushed  density  will  decrease  somewhat  as  the  pressure  behind  the  shock  falls  off. 
Such  behavior  would  cause  the  plot  In  Pig.  6  to  deviate  from  the  straight  line 
character  shown.  Pigs.  7,  8,  9  and  10  are  plots  similar  to  Pigs.  3,  4,  5  and  6 
respectively  for  a  different  set  of  Initial  data. 

The  model  which  we  have  just  described  In  no  way  takes  Into  account  the 
elastic  behavior  of  the  material,  but  assumes  that  the  material  Is  permanently 
crushed.  In  the  real  case,  one  would  expect  the  non-porous  constituent  of  the 
material  to  have  some  elastic  properties  under  compression,  and  while  one  would 
expect  seme  permanent  change  In  the  density  due  to  crushing,  there  might  also  be 
a  responsive  change  In  density  as  a  result  of  compressed  entrapped  air  and  the  e- 
lastlc  properties  of  the  other  constituent  materials.  However,  It  Is  believed 
that  the  large  changes  In  density  will  be  a  result  of  a  crushing  or  compacting 
of  the  material,  particularly  In  the  case  of  porous  media.  As  the  distance  from 
the  source  of  the  explosion  Increases,  the  material  may  not  undergo  crushing  at 
all  but  simply  undergo  an  almost  elastic  expression  to  sxe  new  density  and  then 
upon  release  of  the  applied  forces  be  restored  to  Its  Initial  density.  This  is 
the  behavior  In  the  case  of  seismic  wave  propagation.  In  the  following  section 
we  will  discuss  the  propagation  of  such  elastic  waves. 
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Fig.  7— Position  verbus  tiiae  for  each  of  five  hypothetical  gages  where  the  initial 
configuration  euid  s  j11  behavior  Is  assumed  to  be  given  by  the  following  data. 

Pq  “  3"5  slugs/ft.^;  p/p^  a  2.0;  =  20ft.;  7  a  1.2;  =  10^\b./ft. 

Curves  A,  B,  C,  D  and  E  represent  the  paths  of  the  five  hypothetical  displacement 
gages. 
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Fig.  10— Peak  particle  velocity  measured  at  the  hypothetical 

gages  as  a  function  of  the  initial  distance  of  the  gages 
from  the  explosive  source.  Initial  conditions  are 

p  ■  3. 5  slugs/ft.  p/p  »  2. 0;  a  ■  20ft. ;  7  ■  1. 2; 


10^°lh./ft.^, 
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In  tb*  prtvlouf  snotlon  w«  oonaldtrnd  a  wava  propagating  In  a  soil  which 
pamanantly  oruahad  or  oooqpaotad  to  sona  naw  danalty  different  from  the  Initial 
density.  In  this  saotlon  wa  shall  Investigate  the  opposite  extreme — the  case  where 
the  density  returns  to  Its  Initial  value  after  the  passing  of  the  disturbance. 

The  true  behavior  would  be  expected  to  lie  soiMwhere  between  these  two  extremes 
and  we  will  make  some  attenpt  to  determine  the  type  of  wave  shape  which  we  would 
have  If  the  soil  underwent  a  permanent  deformation  and  also  responded  elastically. 
We  are  now  concerned  with  relatively  small  deformations  as  compared  to  those  con¬ 
sidered  In  the  equations  for  \Mve  propagation  In  crushable  material  and  therefore 
the  equations  describing  seismic  wave  propagation  will  be  linear  and  hence  will 
adsklt  of  closed  form  solutions. 

We  consider  a  homogeneous,  Isotropic  elastic  media  In  which  there  exists 
a  spherical  boundary  between  the  region  of  crushable  behavior  and  that  of  seismic 
behavior.  The  soil  will  be  assumed  to  extend  Infinitely  In  the  radial  direction 
away  from  the  explosive  source.  We  further  assume  that  the  elastic  properties  of 
the  material  are  governed  by  the  usual  Hooke's  law  equations  which  In  the  case  of 
spherical  symmetry  are  given  by 


Oj.  =  (\  +  2p)  (tj.  +  2*^)  -  4^6^' 

e 

y 

(13) 

Og  =  +  2,g)  +  2peg 

e 

(U) 

where  o  is  the  stress  In  the  radial  direction: 
r  ' 


Og  Is  the  stress  In  the  tangential  direction; 

c.  Is  the  strain  In  the  radial  direction; 
r  ’ 


Is  the  strain  In  the  tangential  direction; 
\  and  p  are  the  Lame  parameters. 
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ClASfleal  eqiiwtions  of  motion  in  the  spherically  symmetric  case  can  be 
obtained  in  the  coordinate  system  of  interest  in  a  direct  manner  by  considering 
a  small  conical  voIvom  element  described  as  follovs:  (1)  A  surface  generated 
by  a  solid  angle  from  a  point  in  three  space;  (2)  two  concentric  surface 
generated  by  radii  with  the  aforementioned  point  as  center.  We  assume  the 
following  symmetrical  behavior: 

1.  The  rate  of  change  of  momentum  of  the  conical  volume  Is  equal  to 
the  total  external  force  acting  on  the  volume; 

2.  The  mass  of  the  moving  volume  remains  unchanged  under  deformation 
of  the  medium. 

In  addition,  we  assume  that  the  forces  due  to  the  acceleration  of  gravity  can 
be  neglected  and  that  strains  will  be  relatively  small. 

By  further  consideration  of  the  elemental  volume  previously  described, 
we  c€m  determine  the  relationship  between  the  strain  amd  the  deformation.  These 
relations  are 


S  *  "sr 

,  (15) 

U_ 

«  - — 

9  r 

,  (16) 

where  u.  is  the  radial  deformation: 
r  ' 

r  is  the  radial  distance. 

In  the  classical  fashion,  we  obtain  the  wave  equation, 


C' 


when  \|r  is  a  displacement  or  velocity. 


(17) 
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Where  c  Ib  the  velocity  of  sound  in  the  medium  and  is  given  in  terms  of  the 
l0at  constants  and  the  density.^ 


Ih  the  case  of  small  strains  we  can  write  the  expression  for  the  dilation  as 
a  linear  function  of  the  strain, 

A  -  Cg  +  2€g  .  (18) 

and  from  Eq.s.  (1?)  and  (16)  we  see  that  rA  satisfies  the  wave  equation.  A 
solution  to  the  wave  equation  is  given  by  the  arbitrary  function  ftt  -  (r  -  b^)/c]. 
Eqiuating  this  function  to  rA  we  have 

rA  -  fit  -  (r  -  b^)/c]  ,  (19) 

/ 

where  b^  is  the  radial  dlstsmce  to  the  initial  elastic  boundary. 

We  may  now  write  an  equation  which  includes  the  deformation  in  the  radial 
direction  by  substituting  for  the  terms  and  in  Eq.  (l8)  their  equivalent  forms 
given  in  Eqs.  (15)  and  (l6).  After  some  manipulation  we  obtain  the  partial  differ¬ 
ential  equation, 

(r%)  -  rf  [t  -  (r  -  b^)/c]  =  0  .  (20) 

For  convenience  we  Introduce  a  new  function,  F,  given  by  the  defining  equation, 

^  F[t  -  (r  -  b  )/c]  =  ftt  -  (r  -  b  )/c]  .  (21) 


^TJhomas  T.  Y. ,  Plastic  Flow  and  Fracture  in  Solids,  Academic  Press  (I961) 
p.  60  ff.  Thcmas  also  glve:8  a  more  exact  expression 


where  p  is  the  density,  a  function  of  t  and  r. 
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Ve  than  fom  the  Integral  of  (20)  Incorporating  the  function,  F. defined  In 
Bq,.  (21)  and  after  Integration  by  parte  obtain  the  differential  equation, 

-r^Uy  ■  -  (r  -  b^)/o]  +  cS’[t  -  (r  -  b^)/o]  .  (22) 

Now,  provided  we  know  the  boundary  conditions  for  u^,  that  Is,  if  we  know  the 
defomatlon  as  a  function  of  the  time  and  the  position  of  elastic  boundary, 
we  can  determine  the  deformations  for  all  time  euad  radii.  Similarly,  we  can 
obtain  the  particle  velocity  and  particle  acceleration  for  clII  time, t, at  any 
radius, r,by  differentiation  of  the  expression  for  u^  obtainable  from  Eq.  (22). 

In  addition,  substituting  Eqs.  (1^)  and  (l6)  Into  Eq.  (13)  and  using  the 
relation  Just  obtained  for  u^  In  Eq.  (22),  we  can  determine  for  all  time  and 
radial  distances,  the  radial  stress  wave,  o^,  which  Is  propagating  through 
the  medium: 


r^Cy  -  pcV  +  Vcr  .  (2 

This  eqjiation  ce«i  also  be  used  to  determine  the  f\xnctlon,F,  provided  we  know 
the  pressure  versus  time  on  the  seismic  boundary  or  else  the  pressure  versus 
time  at  some  other  boundary. 

Actually  the  boundary  conditions  required  for  Eqs.  (22)  and  (23)  may 
be  obtained  from  the  results  of  the  computation  carried  out  for  wave  propagation 
in  a  crushable  media  where  we  view  the  sequence  of  events  as  follows.  The 
extremely  high  pressures  generated  within  the  bomb  room  at  detonation  time 
cause  crushing  of  the  soil  out  to  some  radial  distance.  Beyond  this  radius 
the  waves  propagate  seismlcally.  We  may  describe  the  crushed  material  as  a 
large  spherical  piston  which  applies  a  force  to  the  elastically  responding  soil 
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on  ita  extarlor.  This  spherical  piston  laoves  out  to  some  naocimum  radius,  all 
the  vhile  elastically  defoming  the  soil  on  its  exterior. 

At  this  point  In  the  discussion  ve  must  make  some  additional  assumptions 
about  the  structural  Integrity  of  the  thick  vailed  spherical  shell  which  makes 
up  our  piston.  In  general,  we  would  expect  that  the  structural  Integrity  of 
this  piston  would  be  such  that  the  elastic  response  of  the  surrounding 
soil  would  tend  to  drive  the  spherical  piston  back  towards  the  center  of  the 
cavity  until  the  force  due  to  the  cavity  pressure  and  the  restraining  force 
due  to  the  residual  structural  Integrity  of  the  spherical  shell  becomes  suffi¬ 
cient  to  balance  the  static  stress  applied  by  the  surrounding  soil.  Perhaps 
one  might  also  expect  mechanical  oscillation  of  this  spherical  piston  as  a 
conseqiuence  of  the  elastic  properties  of  the  solid  constituent  of  the  crushed 
material. 

For  purposes  of  computation  we  assume  the  following:  the  pressure 
exerted  by  the  spherical  piston  on  the  surrounding  elastic  material  Is  subject 
to  a  step  rise  after  which  the  pressure  relaxes  gradually.  This  would  he 
similar  to  the  effect  we  would  have  if  our  spherical  piston  deformed  the 
surrounding  material  to  some  maximum  value  and  then  partially  releuced, 
maintaining  a  qviasl-static  stress  level  and  therefore  a  permanent  deformation. 

In  order  to  determine  some  qualitative  wave  shapes,  a  solution  was  obtained 
for  Eq.  (23).  The  solution  function  Is  given  by  Eq.  (23)  where  the  function 
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F  is®  T 

F(t)  .  e"*^(A  sin  bT  +  B  cos  bT)  +  |  ^  f  (T)e®^'^  "  '^^/®sln  b(T  -  T)dT  j  (24) 

0 

where  a  *  B^/b^pc 

b  .  2ti(X  +M.)]^^^/bQpc 
T  -  t  -  (r  -  bQ)/o 

f(T)  »  bQP(T,b^)/pc® 

A  and  B  are  constants  determined  by  the  initial  data.  Particular  solutions 
of  Eq.  (23)  have  previously  been  obtained  by  Selberg^,  Goldsmith  and  Allan^,  and 
others. 

We  note  from  the  particle  velocity  plots  of  the  more  distant  hypothetical 
velocity  gages,  which  were  considered  in  Section  II,  that  the  character  of  the 
wave  cem  be  approximated  by  an  exponential  decay  with  a  long  time  constant.  For 


®Kam':e,  E. ,  Differentialgleicbungen  Losungsmethoden  Und  Losungen, 
Chelsea  Publishing  Company,  New  York7  1959»  PP-  412,  413. 

•j 

Selberg,  Henrik  L. .^J'Transient  Compression  Waves  from  Spherical  and 
Cylindrical  Cavities,"  Arklv  For  Fysik,  Vol.  5>  Stockholm,  19$2,  p.  97  ff- 

4 

Goldsmith,  W.  and  Allan,  W.  A. ,  "Graphical  Representation  of  the 
Spherical-Propagation  of  Explosive  Pulses  in  Elastic  Media, "  Journal  of  Acous¬ 
tical  Society  of  America,  Vol.  27,  pp.  47-55,  1955. 


iNOINIIRINO-PHYSICS  CO 


-21  - 


aaludc  dl8turb*no«a  we  have  approxlaately 

P  =  pev  ;  (25) 

where  P  s  pressure; 

p  =  the  density  of  the  soil; 

0  =  propagation  velocity; 

V  =  particle  velocity. 

and  this  Is^lles  that  the  pressure  on  the  boundary  between  the  permanently  crushed 
and  seismic  regions  can  be  approximated  by 

P  =  ^  (26) 

where  t  Is  determined  using  Eqs.  (25)  and  (26)  and  the  data  obtained  in  Sec.  II 
for  the  particle  velocity; 

P^  Is  the  peak  pressure  and  Is  determined  from  Eq.  (25)  and  the  data  of 
Sec.  II. 

We  now  set 

P  =  -Oj.(b^,t  )  ;  (27) 

where  Oy(b^,t  )  Is  the  stress  In  the  radial  direction  at  the  elastic  boundary. 
Substitution  of  this  into  Eq.  (24)  determines  the  function,  F,  which  upon  substi¬ 
tution  into  Eq.  (22)  yields  an  expression  for  the  deformation.  Subsequent  differ¬ 
entiations  yield  expressions  for  the  particle  velocity  and  acceleration. 

Computations  were  ceurrled  out  to  determine  displacement,  velocity,  and 
acceleration  wave  forms.  Two  sets  of  these  computations  are  shown.  One  set  was 
obtained  using  the  elastic  constants  measured  for  salt,  and  the  other  used  average 
values  measured  for  tuff.  The  required  Input  data  Is: 

1.  Young's  modulus,  E; 

2.  Poisson's  ratlo^  u; 
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3.  Dooslt/t  p; 

4«  F*«k  boundary  partiola  veloolty; 

5.  Tima  oonatanty  t; 

6.  Range  to  hypothatioal  gage,  r; 

7.  Radial  distance  to  the  selsnlc  boundary,  b^. 

We  note  that  at  the  wave  front  the  peak  particle  velocity  falls  off  as  1/r.  Sim¬ 
ilarly,  the  defonatlon  falls  off  as  l/r  In  the  case  of  dynamic  loading  but  as  l/r^ 
In  the  case  of  static  loading.  The  displacement  rises  rapidly  to  a  peak  value, 
then  recedes  and  tends  to  oscillate  about  some  steady  state  value  but  damps  out 
rapidly. 

Figs.  11,  12,  and  13  show  the  displacement,  velocity  and  acceleration  as  a 
function  of  time  at  a  hypothetical  gage  300  ft.  from  the  explosive  source.  The 
constants  used  to  describe  the  mechanical  behavior  of  the  soil  are  those  measured 
for  ttiff.  Figs.  14,  15)  and  16  show  similar  curves  for  the  dlsturbcuice  as  measured 
at  a  hypothetical  gage  9OO  ft.  from  the  eiqploslve  source.  Figs.  17  through  22  are 
a  set  of  figures  similar  to  11  through  16,  but  Indicate  the  wave  shapes  where  the 
propagation  medium  Is  sadt.  A  comparison  of  these  two  sets  of  wave  forms  Indicates 
that  those  for  tuff  contain  lower  frequency  coiq>onents  than  the  similar  ones  for 
salt.  Also,  the  displacements  in  salt  are  smaller  thsua  In  tuff.  These  facts  par¬ 
tially  e]q>l«d.n  the  reason  for  shots  in  tuff  being  more  easily  detected  at  a  dis¬ 
tant  station  than  the  same  yield  shot  In  salt. 

Fig.  23  shows  the  dlsplacmnent  wave  measured  at  the  298  meter  station  In 
the  GNOME  event.  The  character  of  the  wave  Is  similar  to  that  predicted  analyti¬ 
cally. 

Fig.  24  shows  the  shape  of  elastic  displacement  waves  generated  by  mechan- 
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Pig.  2l»~-Badlal  distance  of  500  feet  from  the  e3q)loslve  source.  Plots  are  shovn 
for  several  different  values  of  the  time  constant,  t. 
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loftl  disturbUMs  having  various  tias  oonstants.  For  a  oontalnsd  a]q>lo8lon,  tha 
tias  oonstants  art  gaasrally  long  as  indloatad  by  tha  analysis  of  Section  II;  how¬ 
ever^  in  the  oase  of  early  and  rapid  venting  of  the  cavity,  the  time  constant  may 
be  rather  short. 
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IV.  EUSTIC  PRECURSOR 

5 

A  Bo-oalled  elastic  precursor  has  been  observed  In  laboratory  scale  wave 
propagation  experiments  using  naturally  occurlng  earth  materials.  Additionally, 
there  are  some  large  scale  experimental  records^  which  indicate  to  s<mie  extent  the 
existence  of  a  precursor.  Indeed,  one  would  expect  a  precursor  on  theoretical 

n 

grounds  ,  since  experimental  evidence  shows  that  the  shock  wave,  while  propagating 
through  soil  near  the  source  of  the  explosion,  can  vary  in  velocity  from  supersonic 
to  subsonic  and  back  to  sonic.  During  the  time  the  shock  front  is  propagating  sub- 
sonlcally,  we  would  expect  a  sonic  disturbance  to  propagate  ahead  with  an  aii^>litude 
determined  by  the  maximum  elastic  energy  which  can  be  contained  in  the  earth  mater¬ 
ial  under  consideration.  Actually,  the  amount  of  elastic  energy  which  can  be 
stored  in  a  given  mass  of  earth  material  is  subject  to  change  dependent  on  the 
length  of  time  the  energy  is  stored.  Hence,  the  magnitude  of  the  disturbance  which 
is  propagated  through  the  soil  depends  on  the  soil's  stress-strain  history.  For 
the  moment,  we  will  neglect  this  latter  fact  and  set  up  some  simple  relations  to 
describe  the  precursor  in  terms  of: 


^Chabai,  A.  J.,  "Research  and  Analysis  of  Close-In  Data,"  Contract  No. 
DASA-EO-300-60,  Semiannual  Report,  May  1,  1961  -  October  1,  1962,  Sandia  Corpora¬ 
tion  Project  No.  14.041.00. 

^eart,  W.  D.,  "Particle  Motion  Near  A  Nuclear  Detonation  in  Halite,"  Pro¬ 
ject  GNOME,  PNE-108,  Sandia  Cox*poratlon,  Preliminary  Report,  Projects  44.1  and  1.1, 
Feb.  1962. 

7 

Nuckolls,  John  H.,  "A  Computer  Calculation  of  Rainier,"  Lawrence  Radiation 
Laboratory,  Livermore,  PLOWSHARE  SERIES,  Part  I,  Phenomenology  of  Underground  Nu¬ 
clear  Explosions,  UCRL-5675,  May  13-15,  1959. 
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1.  The  velocity  history  of  the  shock  front; 

2.  An  assumed  peak  pressure,  for  elastic  vaves  in  the  medium 
und^  consideration. 

We  take  the  shock  front  velocity  to  he  R  and  the  position  of  the  shock  front 
at  time,  tj  to  be  R. 

In  the  acoustic  approximation,  spherical  divergence  prescribes  that  the 
peak  pressure  be  inversely  proportional  to  r,  the  radial  distance  from  the  point 
of  emission.  We  assume  that  the  pressure,  P^,  is  the  maximum  allowable  pressure 
which  can  be  propagated  acoustically  in  the  medium.  Let  t  be  any  time  during 
which  an  acoustic  disturbance  is  emitted  from  the  traveling  shock  front.  A 
typical  plot  of  time  versus  shock  front  velocity  would  have  the  form  shown  in 
Fig.  25.  Then  the  position,  r^,  of  the  precursor  wave  front  at  some  time,  t, 
where  t  >  t,  is  given  by 


r^  »  R(t)  +  (t  -  t)c  ^  (28) 

where  c  is  the  sonic  velocity  in  the  medium. 

Therefore,  at  this  time,  t,  the  pressure,  P,  of  the  emitted  disturbance  is  given  by 
(in  the  acoustic  approximation) 


P 

Using  the  acoustic  approximation 


PoR(t) 


(29) 


P  »  pcv  ,  (30) 

where  p  is  the  density  of  the  medium 

V  is  the  particle  velocity  in  the  medium 
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Vlg.  2^-k  typical  shock  velocity  versus  time  plot  shovlitg 
supersonic,  subsonic  and  sonic  phases  of  shock  wave 
propagation. 
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we  have 

where  ■  Pg/pc 

Combining  this  result  with  Eq.  (28)  yields 


(31) 


v^R(t) 

""  ‘"kTrl 

Now^  at  a  particular  gage  located  at  a  distance,  r^ 
and  at  a  time,  t,  we  have  from  Eq.  (28) 


(32) 


from  the  explosive  source. 


(33) 


Knowing  the  position  of  the  shock  front  as  a  function  of  time  from  an  analysis 
such  as  that  given  in  Section  II  or  from  experimental  data  we  can  vary  the  parameter, 
T,  and  determine  from  Eqs.  (32)  and  (33) >  the  particle  velocity  as  a  function  of 
time. 

For  t  >  T  the  particle  velocity  history  of  the  existing  wave  is  as 
shown  in  Fig.  26. 

The  pressure  at  the  leading  edge  of  the  elastic  precursor  is  smaller  than 
that  at  the  shock  front  due  to  the  l/r  fall  off  assumed  in  the  acoustic  approximation. 

However,  there  is  yet  another  phenomenon  which  affects  the  wave  form; 
namely,  the  relaxation  of  the  medium  after  passage  of  the  precursor  wave  front. 

This  is  a  property  of  the  medium  under  consideration,  but,  in  general,  the  stress 
is  relieved  internally  in  the  medium  and  thus  the  particle  velocity  decreases  behind 
the  precursor  wave  front  as  illustrated  in  Fig.  27.  A  somewhat  similar  behavior 
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Fig.  26»«Preeursor  1*  tndloatad  bjr  solid  curve.  Nsia  sboek  suv* 
is  shown  by  brohsn  curve.  Curves  ere  not  neeesserllv  to 
scale  since  prscursor  eejr  be  o(  short  duzetioo  eomred 
to  main  distuitance. 


Fig.  27— Qm  effect  of  the  velaaatloa  phsBoswooo  on  the  initial 
precursor  wave  f  oi«. 
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can  be  noted  in  the  case  of  elastic  vaves  previously  described  In  Section  III. 

The  distance  between  the  precursor  wave  front  and  the  shock  front  lengthens 
when  R  <  c,  but  when  the  shock  front  velocity  Increases  to  sonic  velocity,  the 
distance  between  the  fronts  remains  fixed,  with  the  shock  particle  velocity 
dropping  to  v^,  emd  hence  not  distinct  from  the  rest  of  the  wave  form.  Therefore, 
the  final  wave  form  that  propagates  elastically  into  the  medium  has  the  form  given 
by  Fig.  28. 

In  the  event  that  the  shock  front  velocity  drops  only  slightly  below  the 
acoustic  velocity  of  the  soil  or  in  the  event  the  shock  front  velocity  is  sub¬ 
sonic  for  only  a  short  time  the  precursor  may  not  be  detected  by  gages  that  are 
used  primarily  to  detect  the  entire  wave.  Many  field  records  give  no  indication 
of  a  precursor  suggesting  that  the  gage  did  not  respond  to  the  precursor  or  else 
that  the  precursor  in  some  cases  simply  does  not  exist. 

Gages  having  the  proper  response  characteristics  might  be  used  to  obtain 
information  about  the  soil  behavior  near  the  source  of  an  explosion.  For  example, 
if  we  know  the  particle  velocity  as  a  function  of  time  at  a  gage  a  known  distance 
from  the  explosive  source,  as  well  as  the  density,  velocity  of  propagation,  and  the 
maximum  pressure  of  elastic  waves,  then  one  can  determine  from  Eqs.  (30),  (31)  and 
(33)  bhe  new  equations 

pcv(t)r 

R(t) - p - S  ,  (34) 

0 

r  -  R(t) 

T  =  t  -  -  ,  (35) 

which  yield  the  required  shock  wave  velocity  as  a  function  of  time. 
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elastically.  Precursor  may  be  of  considerably  shorter  du: 
compared  to  the  direct  disturbance  than  Is  Indicated  here 
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If  we  had  records  from  several  gages  they  could  he  used  to  obtain  a 
best  fit  for  the  shock  vave  position  and  velocity  as  afvmction  of  time  or  else 
for  each  additional  gage  ve  could  determine  one  of  the  before  required  Inputs 
such  as  density  or  propagation  velocity  of  the  soil. 

Furthermore  with  data  describing  the  position  of  the  shock  front  and 
its  velocity  as  a  function  of  time  ve  can  use  an  analysis  somewhat  the  reverse 
of  that  given  in  Section  II  and  thereby  determine  important  information  about 
the  equation  of  state  of  the  soil  near  the  explosive  source. 


( 
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V.  SPALLING  AND  SBCONDARY  WAVES 

Spalling  of  the  earth  may  result  as  a  secondary  effect  of  an  underground 
nuclear  explosion.  The  following  mechanism  is  the  cause  of  spalling;  the  initial 
shock  wave  generated  hy  the  explosive  source  propagates  towards  the  earth' s  surface 
and  upon  interaction  with  the  surface,  generates  a  rarefaction  wave  which  begins  to 
propagate  back  into  the  soil, having  the  effects  of  increasing  the  particle  velocity 
in  the  direction  normal  to  the  surface  and  at  the  same  time  relieving  the  compressive 
stress  of  the  tail  of  the  direct  pressure  wave.  At  some  depth  the  sum  of;  (l) 
the  stress  of  the  tail  of  the  direct  wave;  (2)  the  lithostatlc  stress;  and,  (3) 
the  tensile  stress  of  the  rarefaction  wave,  may  represent  a  tensile  force  in  excess 
of  that  required  for  rupture  of  the  soil.  At  this  instant  and  at  this  depth,  rupture 
will  occur  and  the  slab  defined  by  the  fractured  surface  and  the  earth's  surface 
will  be  hurled  upward.  There  has  been  ample  experimental  evidence  of  surface  be¬ 
havior  of  this  sort. 

Spalling  has  evidenced  itself  in  the  records  of  vertical  accelerometers 
placed  in  the  soil  above  underground  nuclear  explosions.  For  example,  in  the 
Rainier  shot,  an  acceleration  of  one  gravity  was  measured  by  near  surface  gages  for 
a  period  of  about  l/2  second,  indicating  that  soil  was  being  acted  upon  essentially 
by  gravity  alone  during  that  period.  Such  behavior  was  noted  on  gages  buried  at 
depths  of  200  to  250  feet  below  the  surface,  suggesting  that  rupture  took  place 
beneath  this  depth. 

From  the  VELA-UNIFORM  point  of  view,  secondary  wave  generation  which  results 
from  spalling  can  be  of  considerable  importance,  since  a  large  amount  of  propagated 
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energy  may  be  trapped  in  the  massive  slab  of  spalled  earth  which  eventually  impacts 
the  underlying  earth,  generating  the  observed  secondary  waves.  Depending  upon 
the  geometric  configuration  of  the  spalled  slab  and  the  underlying  surface  it 
impacts,  this  energy  may  represent  a  large  fraction  of  the  hydrodynamic 
energy  produced  by  an  underground  nuclear  explosion.  The  waves  generated  in 
this  fashion  contain  valuable  information  about  the  close-in  soil  characteristics, 
the  wave  sliape  of  the  direct  wave  close-in,  the  depth  of  burst,  and  the  yield. 

Since  these  waves  follow  the  direct  waves  by  a  considerable  length  of  time, 
perhaps  in  the  order  of  one  second  or  more,  the  two  wave  forras  can  be  distinguished. 
By  correlation  of  these  two  types  of  waves,  it  may  be  possible  to  extract  information 
about  the  mechanism  of  their  generation. 

In  the  literature,  one  can  find  numerous  papers  dealing  with  the  mechanism 
of  spalling;  however,  not  a  great  deal  of  information  exists  which  describes  sec¬ 
ondary  wave  generation.  A  phenomenon  akin  to  spalling  in  the  case  of  underwater 
explosions  is  called  bulk  cavitation.  In  this  case,  there  has  been  some  investi¬ 
gation  of  the  generation  of  secondary  waves,  in  addition  to  some  experiments  to 

g 

study  the  phenomenon.  EPCO  is  presently  carrying  out  a  theoretical  investigation 

9 

of  a  three-dimensional  bulk  cavitation  problem  for  the  Office  of  Naval  Research. 
Because  of  the  similarity  between  the  mechanism  of  spalling  and  that  of  bulk  cavi¬ 
tation,  we  find  that  we  may  use  similar  mathematical  models  to  explain  qualitatively 


Tests  conducted  in  June,  I962  by  Underwater  Explosions  Research 
Division  of  the  David  Taylor  Model  Basin.  ! 

^Contract  NONr-3709(00) 
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the  type  of  eecondaxy  waves  ^Ich  are  emitted  in  either  case.  Some  of  the  basic 
differences  which  are  evident  when  the  medium  under  consideration  Is  soil  rather 
than  water,  are  enumerated  below: 

1.  ISxlstence  of  shearing  forces; 

8.  Non-homegenelty; 

3.  Relatively  large  compressibility. 

In  the  latter  portion  of  this  section,  ve  will  undertake  a  rudimentary  analysis 
of  secondary  wave  generation  in  soil. 

For  purposes  of  our  simplified  mathematical  model,  we  will  make  the 
following  assumptions: 

1.  The  earth's  surface  is  completely  flat  and  horizontal  in  the  neigh¬ 
borhood  of  the  explosive  source; 

2.  The  soil  is  homogeneous  and  isotropic; 

3.  The  propagation  velocities  of  both  compression  and  tension  waves 
are  identical; 

4.  The  primary  wave  shape  has  a  sharp  rise  and  an  exponential  talloff ; 

5.  The  maximum  tensile  strength  of  the  soil  is  a  known  constant; 

6.  The  amplitude  of  the  shear  waves  are  negligibly  small  compared 
to  that  of  the  direct  compression. 

While  the  conditions  just  assumed  are  rather  restrictive,  our  model  should  des¬ 
cribe  qjjalitatively  the  features  of  the  secondary  waves  actually  generated. 

It  is  instructive  at  the  outset  to  consider  the  tension  field  which  exists 
when  spallation  does  not  occur.  By  considering  the  sum  of  the  direct  stress  wave, 
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llthostatlc  pressure,  end  the  tensile  stress  caused  by  the  rarefaction  vave  In  a 
three-dimensional  situation  vhere  the  explosive  source  Is  considered  a  point  source, 
ve  obtain  various  typical  stress  fields  where  the  tensile  stress  Is  equal  to  or  in 
excess  of  the  maximum  tensile  strength  of  the  soil.  The  shape  of  field  is  dependent 
upon  the  wave  shape,  depth  of  burst,  yield  and  maximum  tensile  stress;  Typical 
stress  fields  are  shown  In  Fig.  29  for  various  depths  of  burial.  In  the  case  of 
soil,  one  ml{^t  expect  parting  along  some  natural  layer  where  the  tensile  strength 
of  the  soil  Is  less  than  the  magnitude  Indicated. 

For  the  purposes  of  the  qualitative  model  we  assume  that  separation  of  the 
earth  caused  by  spalling  takes  place  along  some  plane  parallel  to  the  earth's 
surface  and,  furthermore,  that  the  spalled  slab  Is  thin  compared  to  the  depth  of 
burst.  At  the  Instant  of  detonation  a  spherical  shock  wave  propagates  radially 
from  the  explosive  source.  For  our  present  purposes,  we  assume  that  the  magni¬ 
tude  of  the  peak  particle  velocity  diminishes  as  the  inverse  first  power  of  the 
distance  from  the  source  of  the  explosion.  Hence  we  have 

u  -  K/R  ,  (36) 

where  R  is  the  radial  distance  from  the  source  of  the  e:qploslon; 
u  is  the  particle  velocity  at  the  dlstsuice  R; 

K  Is  a  proportionality  constcmt. 

The  time,  t  ,  required  for  the  shock  wave  to  reach  the  surface  at  a  distance,  R, 
from  the  explosive  source  Is  given  by 

-  R/c  .  (37) 

where  c  Is  the  propagation  velocity. 

Upon  arrival  of  the  direct  disturbance  at  the  surface,  a  rarefaction  wave  Is 
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surface 


Fig<  29—Region  In  whldi  epalllng  oey  begin  for  the  seme  yield  but 
three  different  depth*  of  bunt.  The  top  figure  Is  for  * 
de^  shot  end  the  bottoa  figure  Is  for  a  staallov  shot.  The 
existence  of  the  ga*  globe  has  been  neglected.  The  contour 
indicates  the  region  under  tensions  In  excess  of  sane  pn> 
determined  value. 
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^  g«n*rated  and  begins  to  propagate  back  Into  the  eoll.  The  rarefaction  wave  having 
a  Btreee  oppoaite  in  sign  to  the  direct  wave  tends  to  relieve  the  stress  in  the 
soil.  In  this  prooesa,  a  portion  of  the  potential  energy  is  converted  to  kinetic 
energy  which  appears  in  the  form  of  increased  particle  velocity.  The  velocity  of 
the  particles  at  the  surface  will  be  directed  normal  to  the  surface  and  have  em 
amplitude  approximately  twice  the  amplitude  of  the  direct  wave  Just  prior  to  im¬ 
pact  with  the  surface. 

With  the  occurrence  of  spalling  we  assume  that  the  resulting  slab  is  hurled 
upward  at  twice  the  particle  velocity  of  the  direct  disturbance  prior  to  impact; 
and  that  it  is  decelerated  by  the  force  of  gravity  alone.  The  flight  time,  t  , 


of  this  spalled  slug  is  then  given  by 


tg  =  4u/g 


where  u^  is  the  vertical  component  of  particle  velocity; 
g  is  the  acceleration  of  gravity. 

The  total  time,  T,  to  impact  as  measured  from  the  lime  of  detonation  is 

T  =  tg  ^  t^ 

or  T  =  R/c  +  4u^g  .  (39) 

In  teims  of  the  particle  velocity,  u,  we  can  write  sm  expression  for  the  vertical 


component,  u 


where  d  is  the  depth  of  burst. 


u^  =  ud/R 


Now,  from  Eq.  (39)  we  can  determine  the  horizontal  distance  at  which  first 


impact  occurs. 
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where  is  the  horizontal  distance  to  the  circle  of  first  Impact. 

This  distance  will  define  a  circle  of  contact  of  the  spalled  slab  with  the  under¬ 
lying  earth.  (In  the  special  case  where  first  impact  occurs  below  surface  zero 
the  circle  of  first  Impact  will  be  a  point. )  In  general  there  will  be  two  circles 
of  closure  or  impact:  one  traveling  Inward  and  one  traveling  outward  from  the 
circle  of  first  Impact.  The  configuration  of  the  Impacting  slab  shortly  after  the 
time  of  first  Impact  is  shown  in  Pig.  30  Plots  of  secondary  wave  arrival  time 
and  impact  velocity  vei'sus  radial  distance  are  shown  in  Fig.  31-  The  velocity 
at  which  this  closure  or  traveling  circle  of  impact  moves  is  given  by  the 

pop 

derivative  of  the  horizontal  distance,  h,  (where  R  »  h  +  d  )  with  respect  to  time, 
t,  and  can  be  obtained  from  Eq.  (39)-  The  horizontal  velocity,  V,  of  closure  is 
found  to  be 

^  ■  V  -  (1  -  .  (42) 

dt  h  gR-’ 

For  large  horizontal  distance,  hence  large  slant  range  R,  it  is  easy  to  see  that 
this  expression  approaches  the  propagation  velocity,  c,  of  the  medium. 

Further  evaluation  of  Bq.  (42)  shows  that  the  inward  traveling  circle 
of  impact  travels  first  at  velocities  well  above  the  seismic  propagation  velocity, 
c,  in  the  medium  but  slows  down  quite  rapidly  and  perhaps  drops  to  velocities  below 
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Fig.  31— Shows  qualitative  results  of  eopputatlons  to  determine  time 
of  impact  and  impact  velocity. 
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J  tiie  propagation  velocity;  finally,  as  the  Inward  traveling  circle  of  liqpact 

approaches  the  point  directly  under  sxirface  zero,  the  horizontal  (inward)  velocity 
of  the  circle  of  Impact  again  Inoreases  to  large  values. 

On  the  other  hand,  the  circle  cf  impact  traveling  outward  from  the  circle 
of  first  Impact  will  also  start  out  at  a  very  high  velocity,  again  above  propa¬ 
gation  velocity  of  the  medium,  gradually  slow  down  to  the  propagation  velocity, 
c,  with  increasing  distance  from  the  explosion. 

The  impacting  of  the  spall  In  this  fashion  produces  a  moving  source  of 
pressure,  and  we  may  now  draw  some  qualitative  conclusions  as  to  the  type  of 
wave  patterns  which  might  be  expected  from  such  a  moving  source.  Using  a  Huygens 
construction  we  may  determine  the  position  of  the  secondary  waves  soon  after  the 
time  of  first  impact.  Fig.  32a  shows  the  shape  of  the  wave  front  shortly  after 
the  time  of  first  impact.  Near  the  circle  of  first  impact  the  closure  or  travel¬ 
ing  impact  point  proceeds  at  supersonic  velocity,  and  hence  a  bow  shock  will  be 
produced.  This  bow  shock  will  be  at  an  angle,  p,  as  shown  in  Fig.  32a,  where  p 
is  given  by 

p  =  arc  sin  (l/M) 

where  M  is  the  Mach  number  of  the  traveling  circle  of  impact.  Fig.  32b  shows  the 
secondary  wave  pattern  well  after  the  time  when  the  inward  traveling  circle  of 
impact  has  fallen  below  the  propagation  velocity  of  the  medium.  The  dotted  lines 
indicated  in  Fig.  32b  denote  gradual  rise  front^  while  the  solid  lines  indicate 
quick  rise  or  shock  fronts.  These  gradual  rise  wave  fronts  occur  because  the 
inward  traveling  circle  of  impact  has  fallen  below  the  medium's  seismic  propagation 

( 
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Fig. 3aa— A  secondary  vave  pattern  soon  after  the  time  of 
first  Impact. 


Fig.  32 1— Secondary  wve  pattern  at  a  time  vhen  the  inwrd  treveling 
circle  of  impeet  has  fallen  below  sonic  veinelty. 
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velocity.  In  Fig.  32b  we  note  that  the  outward  traveling  circle  of  first  Impact 
has  associated  with  It  a  secondary  shock  front  and  that  following  closely  behind 
Is  a  rarefaction  wave  of  very  nearly  equivalent  magnitude.  From  constructions  of 
this  sort  one  can  speculate  on  the  secondary  wave  forms  that  would  be  recorded  at 
any  hypothetical  gage  In  the  field. 

The  foregoing  analysis  and  discussion  Indicate  that  surface  spallation 
phenomena  will  produce- -as  a  consequence  of  the  traveling  point  of  Impact  or 
reloading — secondary  stress  disturbances  which  are  quite  distinct  from  the  body 
waves  emanating  from  an  underground  explosion.  A  more  complete  picture  of 
these  secondary  waves  now  requires  embellishment  of  the  rudimentary  theory  pre¬ 
sented  here,  together  with  numerical  solution  for  configurations  of  interest. 

In  any  further  analysis  the  theory  should  be  extended  to  include  such  factors 
as:  surface  contour;  geologic  layering;  variable  propagation  velocity,  wave 
shape, an.'l  shearing  forces. 
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^  )  VI.  CONCLUSIONS 

The  object  of  this  Investigation  has  been  to  determine  the  wave  shapes 
of  the  seismic  disturbances  that  are  generated  by  an  underground  nuclear  explo¬ 
sion.  The  investigation  took  into  account  (1)  the  generation  of  the  mechanical 
disturbances,  (2)  the  propagation  of  these  disturbances. 

Two  means  of  wave  generation  have  been  considered:  (a)  generation  of  waves 
from  the  expanding  cavity  (body  waves) ;  (b)  generation  of  waves  due  to  surface 
spall  and  subsequent  Impact.  Additionally  we  have  considered  two  types  of  wave 
propagation:  (a)  waves  in  crushable  material  and  (b)  waves  in  elastic  materlsil. 

These  two  means  of  propagation  might  be  thought  of  as  the  two  extreme  ways  In 
which  waves  cein  propagate  in  solid  materials. 

The  elastic  precursor,  a  wave  shape  detail  of  the  direct  or  body  wave  dis¬ 
turbance,  has  been  simply  treated  so  as  to  show  the  the  possible  gains  to  be  ob¬ 
tained  from  a  more  detailed  study. 

In  the  foregoing  analysis,  we  have  observed  the  following: 

1,  Details  of  the  precursor  may  be  significant  in  determining  the 
behavior  of  the  soil  (details  about  the  equation  of  state)  in  the 
region  close  in  to  a  nuclear  explosion.  Also,  this  disturbance  is 
the  first  received  by  a  seismic  station  and,  hence,  will  not  be 
aubject  to  interference  from  mechanical  disturbances  propagating 
along  other  paths  from  the  same  source.  Therefore,  more  emphasis 
should  be  placed  on  instrumentation  to  determine  the  character  of 
the  leading  edge,  and  on  analysis  to  determine  the  expected  shape 
of  the  leading  edge. 

2.  Waves  generated  by  surface  spall  and  subsequent  impact  of  the  spall 
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may  be  used  to  gain  Inforoation  about  the  environment  of  the  explo¬ 
sive  source  and,  in  particular,  about  the  close-in  geonetry.  Com¬ 
parison  of  these  secondary  waves  with  the  direct  disturbance  may 
yield  important  information  about  the  mechanism  of  wave  generation. 

3.  A  simple  crushable  equation  of  state  can  be  used  in  conjvmotion  with 
equations  of  motion  to  determine  the  gross  features  of  wave  propa¬ 
gation  in  the  crushable  region  near  an  explosive  Source.  A  more  ac¬ 
curate  analysis,  if  needed,  should  include  the  effect  of  variable 
density  behind  the  wave  front.  Shearing  forces  should  also  be  in¬ 
cluded  in  a  more  profound  analysis,  especially  in  respect  to  waves 
propagating  through  the  region  near  the  seismic  transition  zone. 

In  the  very  close-in  region,  the  analysis  may  need  to  include  ther¬ 
modynamic  theory  although,  in  the  crushable  region  of  our  interest, 
it  is  believed  that  such  an  inclusion  would  cause  only  a  perturba¬ 
tion  to  the  resiilts. 

4.  In  the  extension  to  the  work  reported  herein,  the  authors  are 

developing  a  general  theoretical  framework  which  covers  the  elastic 
and  crushable  wave  analysis  as  special  cases— and^  importantly, 
should  permit  investigation  between  these  two  extremes. 
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